Tilt-angle landscapes and temperature dependence of the conductance in 
biphenyl-dithiol single-molecule junctions 
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Using a density-functional-based transport method we study the conduction properties of several 
biphenyl-derived dithiol (BPDDT) molecules wired to gold electrodes. The BPDDT molecules 
differ in their side groups, which control the degree of conjugation of the 7r-electron system. We 
have analyzed the dependence of the low-bias zero-temperature conductance on the tilt angle ip 
between the two phenyl ring units, and find that it follows closely a cos^ ip law, as expected from an 
effective ir-orbital coupling model. We show that the tilting of the phenyl rings results in a decrease 
of the zero-temperature conductance by roughly two orders of magnitude, when going from a planar 
conformation to a configuration in which the rings are perpendicular. In addition we demonstrate 
that the side groups, apart from determining (p, have no influence on the conductance. All this is in 
agreement with the recent experiment by Venkataraman et al. [Nature 442, 904 (2006)]. Finally, we 
study the temperature dependence of both the conductance and its fluctuations and flnd qualitative 
differences between the examined molecules. In this analysis we consider two contributions to the 
temperature behavior, one coming from the Fermi functions and the other one from a thermal 
average over different contact conflgurations. We illustrate that the fluctuations of the conductance 
due to temperature-induced changes in the geometric structure of the molecule can be reduced by 
an appropriate design. 

PACS numbers: 73.63.-b, 73.63.Rt 
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I. INTRODUCTION 

In atomic-scale conductors the precise positions of the 
atoms have a decisive influence on the electronic trans- 
port properties.ii'^'^ In the case of metal-molecule-metal 
contacts the importance of such details often complicates 
the reproducibility of the experimental resultsi^i^ For 
this reason statistical analyses of the experimental data, 
such as conductance histograms, have become indispens- 
able for exploring the charge-transport characteristics of 
single-molecule contactsi^'i'^'^'i^ 

Recently, making use of conductance histograms, 
Venkataraman et al. were able to reveal the influence 
of molecular conformation on the conductance of single- 
molecule contactsjii In their experiments, these authors 
investigated biphenyl-derived molecules, where different 
side groups were used to control the tilt angle ip between 
two phenyl rings. Thereby, the extent of the delocal- 
ized TT-electron system of the molecules could be varied. 
They found that the conductance exhibited a character- 
istic cos^ (fi behavior, as expected from a simple effective 
TT-orbital coupling model i^^i^^ 

Motivated by the experiment of Ref. [ll|, we analyze 
theoretically the charge-transport properties of three dif- 
ferent biphenyl-derived dithiol (BPDDT) molecules con- 
nected to gold electrodes. For simplicity, we refer to 
these molecules as R2, S2, and D2 (Fig. [I|. While R2 
is the conventional biphenyl, the other two molecules, 
2,2'-dimethyl-biphenyl (S2) and 2,6,2',6'-tetramethyl- 




Figure 1: (Color online) Biphenyl molecules R2, S2, and D2. 
For S2 one hydrogen atom in the ortho position with respect 
to the ring-connecting carbons in each phenyl ring of R2 is 
replaced by a methyl group. For D2 also the second ortho- 
positioned hydrogen is substituted by a methyl group. 



biphenyl (D2), differ from R2 by the incorporation of one 
or two methyl groups in the ortho position with respect 
to the ring-connecting carbon atoms. 

As an important difference to Ref. [ll|; we bind the 
molecules R2, S2, and D2 to gold using thiol end groups 
instead of amino groups. It has recently been shown that 
amino groups are better suited to establish a reproducible 
contact of a molecule to a gold electrodei^i^ These find- 
ings have been explained by the less directional character 
of the amine-gold bond as compared to the thiol-gold 
linkage. Nevertheless, thiol groups remain a frequent 
choice to establish the electrode- molecule contact (^li^ii^ 
and the molecules R2, S2, and D2 with acetyl-protected 
sulfur end groups have recently been synthesized.— Be- 
sides, it is the internal structure of the molecules that is 
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most important for the charge-transport characteristics 
discussed below. 

In this work we investigate the effect of the degree of 
TT conjugation on the conductance of BPDDT molecules 
connected to gold electrodes. For this purpose, we de- 
scribe the electronic structure of the single-molecule con- 
tacts at the level of density functional theory (DFT). 
We demonstrate that, in agreement with the experi- 
ments of Ref. [Ill, a cos^ ip behavior of the low-bias zero- 
temperature conductance is obtained. This behavior is 
by and large independent of the methyl side groups intro- 
duced. We find that the breaking of the conjugation re- 
duces the zero-temperature conductance by roughly two 
orders of magnitude. In addition, we study the tempera- 
ture dependence of the conductance for all our molecules. 
For this we take two contributions into account. The first 
one comes from the broadening of the Fermi functions 
of the leads, the other one from a thermal average over 
different geometric configurations. In our analysis we ob- 
serve that S2 and D2 exhibit a monotonously increasing 
conductance as a function of temperature, while for R2 
the temperature dependence is non-monotonous. Finally, 
we demonstrate that the temperature-fiuctuations of the 
conductance of single-molecule contacts can be reduced 
by an appropriate design of the geometric structure of 
a molecule. This design should aim at stabilizing the 
molecule with respect to the internal degrees of freedom 
that are most relevant for its conduction properties. Due 
to the elimination of uncertainties about the molecule's 
internal structure in a contact, a more reliable compari- 
son of experimental and theoretical results on the charge- 
transport characteristics of single-molecule junctions can 
be expected. 

The rest of the paper is organized as follows. In Sec.lTTl 
we outline the methods used to compute the electronic 
structure, geometry, and conductance of the molecular 
contacts discussed below. Sec. Hill is devoted to the dis- 
cussion of the results for the conductance of the three 
BPDDT molecules, in particular their tilt-angle depen- 
dence and temperature behavior. Technical details on 
these issues are deferred to Apps. [A] and [BJ Finally, we 
summarize our results in the conclusions. Sec. IIVI 



II. THEORETICAL MODEL 

In this section we present the methods applied in 
our work. These include the procedures for comput- 
ing the electronic structure, the contact geometries, and 
the conduction properties of the molecular contacts. For 
further details on our method we refer the reader to 

Refs.[ii[i3[il. 



A. Electronic structure and contact geometries 

For the determination of the electronic structure we 
employ DFT as implemented in the RI-DFT module 




Figure 2: (Color online) Molecules R2, S2, and D2 contacted 
at both ends to Au [111] electrodes via a sulfur atom in a 
three-fold binding position. Faintly overlayed on the ground- 
state structure of the single-molecule contacts are geometries, 
where the left ring is rotated to if — ^po ± 30°. Here v?, as 
indicated for R2, is the tilt angle between the planes of the two 
phenyl rings, and (/Jo is its ground-state value. The division 
of the junctions into the left (L), central (C), and right {R) 
regions, of relevance in the conductance calculations, is also 
shown. 



of the quantum chemistry package Turbomole v5.7 
(Refs. [lli^). All our calculations, including the elec- 
trode description, are done within TuRBOMOLE's stan- 
dard Gaussian basis set, which is of split-valence qual- 
ity with polarization functions on all non-hydrogen 
atoms . As the exchange-correlation functional we 
use BP86 (Refs. H^IHIll). All calculations were per- 
formed in a closed-shell formalism, and total energies 
were converged to a precision of better than 10~^ atomic 
units. 

Our contact geometries are displayed in Fig. [21 For 
their determination we calculate at first the geometric 
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structure of a gold (Au) [111] pyramid with a thiolated 
benzene molecule on top. This pyramid consists of three 
layers of Au with 3, 6, and 10 atoms. The tip atom of 
the pyramid is missing so that the sulfur (S) atom of 
the benzene binds to a threefold hollow site on top of 
the Au structure. We relax all atomic positions except 
for the layers containing 6 and 10 atoms, which are kept 
fixed. In particular, their lattice constant is set to 4.08 A, 
the experimental bulk atomic distance of gold. Next, we 
compute the geometry of the biphenyl molecules R2, S2, 
and D2 of Fig. [T] with the hydrogen atoms in the po- 
sition 4 and 4' substituted by sulfur atoms, which are 
bonded to a single gold atom, respectively. For each of 
these molecules we replace the single Au atoms on each 
side with the previously mentioned Au [111] pyramids, 
where the thiolated benzene molecule has been removed. 
In this process we take care that the sulfur atoms of the 
biphenyl molecules are in the old positions of the sulfur 
atoms of the thiolated benzene on top of the Au pyra- 
mids, and that the S — S molecular axis and the [111] 
direction are aligned. The ground-state contact geom- 
etry is subsequently obtained by relaxing the complete 
structure, where we keep only the terminal two gold lay- 
ers on each side fixed (atoms shown in blue in Fig. [2). 
As above, the lattice constant in these layers is 4.08 A. 
The relaxations thus include all atoms of the molecule, 
two sulfur atoms, and six Au atoms, three on each side. 
In the determination of the contact geometries we let the 
maximum norm of the Cartesian gradient decay to values 
below 10~^ atomic units. 



B. Transmission function 

To compute the charge transport we apply a method 
based on standard Green's function techniques and the 
Landauer formula expressed in a local nonorthogonal 
basis.— i^>21 The local basis allows us to partition the ba- 
sis states into left (L), central (C), and right (i?) parts, 
according to a division of the contact geometry. Thus 
the Hamiltonian (or single-particle Fock) matrix H, and 
analogously the overlap matrix S, can be written in the 
block form 





( Hll 









HcL 


Hcc 


HcR 




K 


Hue 





Within the Green's function method the energy- 
dependent transmission t{E) is expressed as2^ 

T(E)^Tr[TLG},c^RGlc]. (1) 
with the Green's functions 

G},c{E) = [EScc - Hcc ~ ^1{E) ~ S?,,(£;)]-i 
and Gqc — [Gqq]\ the self energies 

^x{E) = (Hex - EScx) 9xxiE) [Hxc — ESxc) , 

(2) 



the scattering rate matrices TxiE) = — 2Im [S^(£')], 
and the electrode Green's function gxx^ where X = 
L, R. 

In Fig. [2] we show how we divide our contacts into 
the L, C, and R regions. In this way we obtain Hcc 
and Sec for the C region, which consists of the BPDDT 
molecule and three gold atoms on each side of the junc- 
tion. The L and R regions are made up of the two ter- 
minal gold layers on each side of the gold pyramids (blue 
shaded atoms of Fig. [2]) , which have been kept fixed to 
bulk atomic distances in the geometry relaxations. The 
matrices Hex and Sex, extracted from these finite con- 
tact geometries, serve as the couplings to the electrodes 
in the construction of T,x{E)- However, the electrode 
Green's functions QxxiE) in Eq. ^ are modeled as sur- 
face Green's functions of ideal semi-infinite electrodes. 
In order to obtain these surface Green's functions, we 
compute the electronic structure of a spherical gold clus- 
ter with 429 atoms. From this we extract the Hamilto- 
nian and overlap matrix elements connecting the atom 
in the origin of the cluster to all its neighbors and, using 
these "bulk parameters", construct a semi-infinite crystal 
which is infinitely extended perpendicular to the trans- 
port direction. The surface Green's functions are then 
calculated from this crystal using the so-called decima- 
tion technique.— We have checked that the electrode con- 
struction in the employed nonorthogonal basis set has 
converged with respect to the size of the Au cluster, 
from which we extract our parameters.— In this way 
we describe the whole system consistently within DFT, 
using the same nonorthogonal basis set and exchange- 
correlation functional everywhere. 

We assume the Fermi energy Ep to be fixed by the 
gold leads. From the AU429 cluster we obtain a Fermi 
energy for gold of Ep — —5.0 eV, which is chosen to 
lie halfway between the levels of the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) of the cluster of —4.96 and 
—5.01 eV, respectively. 



C. Conductance 

The low-bias conductance in the Landauer formalism 
is given b y^^i^° 
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In this expression T is the temperature. Go — 
2e^/h is the conductance quantum, t^{E) is the 
energy-dependent transmission [Eq. ([T])], and f{E) — 
1/ \^Q-(E-EF)/kBT _^ ^Yie Fermi function with Boltz- 
mann's constant ks- The factor —df{E)/dE is the so- 
called thermal broadening functioni^^ For zero temper- 
ature, Eq. im reduces to G^{T = OK) = Gqt^{Ef). 
Here and henceforth we index G and r with (/?, which pa- 
rameterizes different geometrical contact configurations. 
In our case these configurations correspond to different 
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molecular conformations with the tilt angle ip between 
two phenyl rings (Fig. [5]). At finite temperature, tilt an- 
gles If differing from the minimum-energy, ground-state 
value can be accessed. We account for this additional 
temperature-dependent contribution to the conductance 
by the thermal average 



.,30.31 



G(T) = (G^(r))^ 



(4) 



with (■•■)^ = /d(^e~^-A-«^ (••■)// (/(pe-^^/'^^^. In 
this expression is the energy of the metal-molecule- 
metal contact for angle Lp. The contribution of G^(T) to 
G{T) can be seen as a "electronic" or "lead-induced" tem- 
perature dependence, because it follows from the Fermi 
functions of the electrodes. Its determination is discussed 
in App. [Al On the other hand, the Lp average represents 
a "configuration- induced" contribution to G{T). For the 
later discussion, we also introduce the variance 



5G{T) = M{G^{T)~G{T)Y)^ 



that describes the fiuctuations of the conductance. 



(5) 
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Figure 3: (Color online) Transmission T^g{E) as a func- 
tion of energy E for the ground-state geometries of the con- 
tacts shown in Fig. [2] The zero-temperature conductances 
G^o(r = OK) of molecules R2, S2, and D2 are 9.2 x IQ-^Go, 
1.9 X 10~"'Go, and 1.2 x lO^^'Go, respectively. The vertical 
dashed line indicates the Fermi energy Ef. 



III. RESULTS AND DISCUSSION 

Let us first discuss some properties of the isolated 
molecules (Fig.[T]). For R2, S2, and D2 we obtain phenyl- 
ring tilt angles p of 36.4°, 90.0°, and 90.0°, respec- 
tively. The tilt angle of R2 is a result of the interplay 
between the tt conjugation, which tries to flatten the 
structure (ip = 0°), and the steric repulsion of the hy- 
drogen atoms in the ortho positions with respect to the 
ring-connecting carbons, which favors tilt angles close 
to p — 90° (Ref. Is^ ). The methyl groups introduced 
in S2 and D2 increase the steric repulsion and cause a 
larger p. As a consequence, the conjugation between the 
phenyl rings is largely broken in S2 and D2, whereas R2 
still preserves a reasonable degree of delocalization of the 
TT-electron system over the whole molecule. This fact is 
clearly reflected in the change of the HOMO-LUMO gaps 
A, which are 3.85 eV for R2, 4.74 eV for S2, and 4.70 eV 
for D2. Thus, A increases by roughly 1 eV when going 
from R2 to S2 or D2. This suggests that the molecules 
S2 and D2 will show a more insulating behavior than R2, 
when they are incorporated into a molecular contact. 

Now, we study the geometric structure of the metal- 
molecule-metal contacts. In Fig. [2] we show the biphenyl 
molecules contacted at both ends to gold electrodes via 
sulfur bonds, where the sulfur resides on the threefold 
hollow position of Au [111] pyramids. The molecular con- 
formation in the junction is very similar to the ground- 
state structure of the isolated molecule, as there is no in- 
ternal stress on the molecule in this binding position.— 
In particular, we obtain ground-state tilt angles ipQ of 
33.8°, 89.3°, and 89.7° for R2, S2, and D2, respectively. 

In order to analyze the conduction properties of these 
molecular junctions, we have computed the transmission 



T^p^{E) as a function of energy for the ground-state ge- 
ometries of the contacts (angle piQ in Fig. The re- 
sults are plotted in Fig.O and our transmission curve for 
R2 agrees well with previous theoretical studies.— Obvi- 
ously, T^g{E) is dominated by a gap around the Fermi 
energy Ep, which reflects the HOMO-LUMO gaps A of 
the isolated molecules. As can be expected due to the 
similar geometric conformations of molecules S2 and D2 
with pq ~ 90°, their transmission curves closely resemble 
each other. However, the most important observation to 
be made from Fig. [3] is the great reduction of the trans- 
mission t^q{Ef) at the Fermi energy for S2 and D2 as 
compared to R2. In particular, the conductance of S2 
(D2) is lower than that of R2 by a factor of 48 (77), 
i.e. roughly by two orders of magnitude. This clearly re- 
veals the importance of the conjugated 7r-electron system 
for the charge transport in biphenyl molecules,^- In ad- 
dition, it shows that the conductance can be tailored by 
means of adequate side groups that force the biphenyl 
molecules to adopt different molecular conformations.— 
To investigate the dependence of the conductance on 
the tilt angle in more detail, we have continuously varied 
p for all the contacts. We do this by rotating one of the 
phenyl rings with respect to the other, as illustrated in 
Fig. [21 without relaxing the contact geometries for tilt 
angles deviating from p^. By changing p, we obtain the 
results depicted in Fig. [H where the total energy E^ and 
the conductance G^ (T = K) are plotted as a function 
of p (Ref. i35 .36 . 37 ) . We explore the tilt-angle intervals 
of [0°, 360°[ few R2, [60°, 300°] for S2, and [60°, 120°] 
for D2 (Ref. [3^ . In each case the angular resolution is 
A(^ = 2°. 

In the energy curve E^ of molecule R2 there are eight 
extrema visible, four minima and four maxima. They 
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Figure 4: (Color online) Landscapes of the total energy E^p (upper panels) and the conductance G^{T = OK) (lower panels) 
as a function of tilt angle (f> for the molecules R2, S2, and D2. Dotted vertical lines indicate positions of extrema in for 
the respective molecules (all panels). In the lower panels a function of the form a + (icos^ tp has been fitted to G^{T = OK) 
(see the legend). For the fit parameters we obtain a = 5.95 x 10~^Go and /3 = 1.35 x 10"^Go (R2), a = 1.81 x 10"''Go and 
P = 1.35 x 10~^Go (S2), and a = 1.58 x 10~*Go and P = 1.42 x 10"^Go (D2). 



are located at 34°, 144°, 214°, and 324° for the energy 
minima and 90°, 178°, 270°, and 358° for the maxima. 
Due to the symmetry of the biphenyl molecule, one would 
expect a 180° periodicity and a mirror symmetry of both 
and G^{T = OK) with respect to 0° (or equivalently 
90°, 180°, or 270°). While the 180° periodicity is present 
for Eip, the mirror symmetry is violated, as one can see 
in Fig. m The reason for this is that the hydrogen atoms 
have been fixed in their positions with respect to the 
phenyl rings as obtained for the ground-state tilt angle 
(fio- They are standing slightly away from the phenyl 
ring planes in this position, which leads to the observed 
violation of the mirror symmetry.— Contrary to E^p , all 
expected symmetries are restored for the conductance. 
In particular, G^{T = OK) possesses only two minima 
at 90° and 270° and two maxima at 0° and 180°. As a 
function of tilt angle, the conductance changes from 2.0 x 
10~''Go in the minima to 1.4 x 10~^G'o in the maxima, 
that is, it changes by a factor of 70. 

In the case of molecules S2 and D2, G^{T = OK) fol- 
lows closely the shape of the energy curve. For S2 there 
are two minima in E^p at 94° and 268° with conductances 
of 2.2 X 10~''G'o, separated by a local maximum at 174° 
with a conductance of 1.3 x 10~^G'o. This corresponds to 
a ratio of 60 between the maximum and minimum con- 
ductance. D2 exhibits an energy minimum at 90° and 
the conductance at this point is 1.2 x 10~^G'o. 



The close agreement of the minimal conductances for 
R2, S2, and D2 (Figs. [3] and [U is remarkable. In the con- 
ductance minima the conformations of these molecules 
are the same, except for their different side groups and 
their slightly varying orientations with respect to the gold 
electrodes. These observations demonstrate that the side 
groups control the conformation, but otherwise have little 
impact on the zero-temperature conductance. This is in 
agreement with the experimental observations of Ref . [ill. 

The large ratios between maximal and minimal con- 
ductances (70 for R2 and 60 for S2) reported above 
highlight the relevance of the extent of the conjugated 
TT-electron system on the conduction properties of the 
biphenyl molecules. In order to further investigate this, 
we have fitted the Gp{T = OK) curves of Fig.|4]to func- 
tions of the form a + {3 cos^ (see the figure caption 
for the obtained fit parameters). A behavior G^{T — 
OK)/G'o oc cos^ (/3 is expected if the coupling between 
the TT-electron systems of the two phenyl rings plays the 
dominant role in charge transport, as discussed in more 
detail in App.[Bj For all three molecules our fit matches 
Gip{T = OK) very well. What is more, we obtain a 
very similar parameter for all of them. On the other 
hand, the small but positive values of a indicate that 
the conductance at perpendicular tilt angles {(p = 90° or 
270°) does not vanish entirely, as a pure cos^ (p depen- 
dence would suggest. This observation was also made in 
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Figure 5: (Color online) Behavior of the conductances G{T) 
and {G^{T = 0K))_^ as a function of temperature T for the 
molecules R2, S2, and D2. 



Ref. luJ. The absence of a complete blockade of the trans- 
port can be understood by the presence of other than the 
TT-TT couplings. 

Next, we analyze the behavior of the conductance with 
respect to temperature. In addition to G(T) = {G^{T))^ 
[Eq. we study (G^(r = OK))^. In this way we 
can quantify the lead-induced contribution to the tem- 
perature dependence of G{T). To perform the average 
{Gip (T = K))j^ we use the energy and conductance land- 
scapes E^p and G^p{T = OK) of the gold-molecule-gold 
contacts as shown in Fig. IH For G{T), instead, we have 
calculated the transmission function for each angle in an 
interval around Ep in order to obtain G^{T) (see the 
explanations in SecHlland App.[A]).— 

The temperature-dependent conductances G(T) and 
{Gy,{T = OK))^ are plotted in Fig. [5] for the molecules 
R2, S2, and D2 for temperatures T between and 300 K. 
With respect to the behavior of G{T) we observe qualita- 
tive differences for the three molecules considered. S2 and 
D2 exhibit a monotonously increasing conductance G{T) 
with increasing T. In contrast, after an initial drop, a 




50 100 150 200 250 300 
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Figure 6: (Color online) Fluctuations SG{T) of the conduc- 
tance as a function of temperature T for the molecules R2, 
S2, and D2. 



non-monotonous weak temperature dependence is found 
for R2. 

The small differences between G{T) and 
{G^{T = OK)) ^ for S2 and D2 indicate that for 
these molecules the lead-induced contribution to the 
temperature dependence can be neglected as compared 
to the configurational one. The monotonous increase 



of G{T) can therefore be understood by E^ and 
G^{T = OK) (Fig. [4]). The ground-state or equivalently 
zero-temperature configurations for both molecules 
correspond to conformations with minimal conduc- 
tances. Therefore elevated temperatures give access to 
conformations with higher conductance values, resulting 
in the observed steady increase of G(T). For molecule 
R2 the situation is different. Here, the energy E^^ at the 
ground-state tilt angle of (/?o = 34° does not correspond 
to a minimum of Gp(T = OK). Elevated temperatures 
give access to both higher and lower conductances and, 
as a result, {G^{T = 0K))_^ exhibits no clear trend. The 
differences between G(T) and {Gip{T — QK))^ signify 
that for R2 both contributions to the temperature 
dependence of the conductance, namely the lead-induced 
and the configuration-induced ones, play an equally 
important role4^ 

Finally we analyze the fiuctuations of the conductance 
5G{T) [Eq. lO], which we have plotted in Fig.[6]for tem- 
peratures T ranging between and 300 K. In each case 
5G{T) increases monotonously with T. This is clear, 
as finite temperatures give access to conductances dif- 
fering from the ground-state conductance G^^ (T = K) 
(Fig. SI). It is also evident from Fig. [6] that 5G{T) is 
largest for R2 and smallest for D2. Indeed, because 
of its two methyl side groups, D2 is the most rigid of 
the three molecules with respect to ring tilts, while R2 
can access a large range of conductance values due to 
its shallow energy landscape (Fig. [5]). Since S2 has 
only one methyl side group, it is an intermediate case. 
From Fig. [6] we obtain the temperature-fiuctuation ratios 
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5GYi2{T)/5Gs2{T) = 4.2 and 5Gs2{T) / 5G^2{T) = 4.2 
for T = 300 K, where 5Gx2{T) refers to the variance of 
the molecule X2. 

In experiments with molecular junctions the conduc- 
tance in the last plateau of an opening curve is gener- 
ally attributed to that of a single molecule. In practice 
the measured conductances are always time averages over 
any fast fluctuations of the contact geometry, in particu- 
lar the internal conformations of the molecule. In terms 
of our definitions, G{T) represents such a time average 
for a given contact realization, while 5G{T) describes the 
fast fluctuations. 

The time-averaged single-molecule conductance may 
vary considerably from one junction realization to 
another.— 1^1^ These variations, and hence the peak widths 
in conductance histograms, can be attributed to two 
types of factors. The flrst one is due to changes at the 
molecule-electrode interface and in the contact environ- 
ment, and the second one is due to modiflcations of the 
molecule's internal geometric structure. Concerning the 
first contribution, the surfaces of the metallic electrodes 
are atomically rough and disordered, and the molecule 
binds differently to the electrodes in every realization of 
the junction.— As a result, the interface-related vari- 
ability of the single-molecule conductance is hard to con- 
trol at present, although more reproducible results can 
be achieved by a proper choice of the binding groups.— 
Regarding the second point, however, the recent pos- 
sibilities of chemical synthesis allow the structure of a 
molecule to be designed. In order to measure the con- 
ductance of a single molecule more reproducibly, such 
a design should aim at making the molecule "rigid" (for 
example by means of side groups, as in the examples con- 
sidered above) i^^^i^ In this way the variability due to 
the changes in the internal structure is reduced, because 
the bonded molecules stay closer to their ground-state 
conformations in isolation. 

In our analysis several simplifying assumptions have 
been made. In particular, we have concentrated on a 
certain realization of a single-molecule junction (Fig. [2]) 
and all temperature-induced changes at the electrode- 
molecule interface have been neglected. Furthermore, 
only one configurational degree of freedom of the 
molecule, the tilt angle ip, has been taken into account, 
and we treated it as a classical variable. Nevertheless, 
our analysis serves to illustrate the importance of the 
temperature-related effects on the average conductance 
and its fluctuations, and how these can be controlled 
by an appropriate design of the molecules. By reduc- 
ing uncertainties about the contact geometries in this 
way, comparisons between theoretical and experimental 
results can be made with a higher degree of confldence. 



IV. CONCLUSIONS 

In conclusion, we studied the charge-transport prop- 
erties of different dithiolated biphenyl derivatives. We 



showed by means of density-functional-based methods 
that the conduction properties of these molecules are 
dominated by the degree of 7r-electron delocalization. A 
broken conjugation, induced by side groups, was found 
to suppress the conductance by roughly two orders of 
magnitude. By varying the tilt angle (p between the dif- 
ferent phenyl rings, we observed a clear cos^ (p behavior 
of the zero-temperature conductance. However, the sup- 
pression of the conductance for perpendicular ring con- 
flgurations is not complete due to the presence of other 
than TT-TT couplings between the phenyl rings. We showed 
that the methyl side groups in the biphenyl molecules 
control the conformation, but they have little impact on 
the zero-temperature conductance otherwise. All these 
flndings are in agreement with the experimental results 
of Ref . mi. 

Based on the energy landscapes with respect to ring 
tilts, we also determined the temperature dependence of 
the conductance. Here we considered two different con- 
tributions. The flrst one originates from the Fermi func- 
tions of the leads, while the other one is due to a ther- 
mal average over different contact configurations. We 
observed qualitatively different temperature characteris- 
tics for the well-conjugated biphenyl molecule as com- 
pared to the molecules whose conjugation was broken 
by means of methyl side groups. Furthermore, we il- 
lustrated that an appropriate design can help to reduce 
temperature-induced conductance fluctuations by stabi- 
lizing a molecule in a conformation close to its ground- 
state structure in isolation. In this way uncertainties 
with respect to the molecule's internal structure are re- 
duced, and a more reliable comparison between theoret- 
ically and experimentally determined charge-transport 
properties of single-molecule junctions seems possible. 
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Appendix A: TEMPERATURE DEPENDENCE OF 
THE CONDUCTANCE 

In order to evaluate the temperature behavior of the 
conductance G^{T) as defined in Eq. we make a 
Sommerfeld expansion.— For this we compute t^{E) 
[Eq. I!])] at 11 equally spaced points in the energy in- 
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terval -0.5 eV < E — Ep < 0.5 eV around the Fermi 
energy Ep for every tilt angle ip. This interval is chosen 
large enough that the thermal broadening function^^ 



b(E,T) = -^f(E) = ^— 



sech 



E — Ep 
2kBT 



has decayed to small values even for the highest tem- 
peratures considered. (For example one gets b{Ep ± 
0.5eV,300K)/6(i;i.,300K) = 1.59 x 10"*.) Next, we fit 
a polynomial of order N = 10 to t^{E) at these energy 
points according to 



N 



With the coefficients 4"^ = cTr^iE) / dE'^\j^^j^^ /n\ de- 
termined from the least squares fit [where in particular 



{E - Epf . 



T^{Ep)], the temperature dependence of G^{T) is 



given as 



G^{T) = Go 



2m 



xC(2to) (2m)!r(2™) (fc^T) 



In this expression \_N/2\ is the largest integer smaller 
than or equal to N/2 and C,{x) is the Riemann zeta func- 
tion. 



Appendix B: EFFECTIVE vr-ORBITAL 
COUPLING MODEL 



The dependence of charge transfer on the tilt angle ip 
between two phenyl rings has been inspected previously 
in Refs. fl^flSl In this appendix we discuss explicitly, how 
the cos^ </? behavior of the conductance can be understood 
based on an effective 7r-orbital coupling model within the 
Green's function formalism. 

For this purpose we bring the transmission function 
Tio(E) [Eg. ^] into a sHghtly different form, following 
Ref. [l3 We assume that the C part of our contacts can 
be divided into two regions 1 and 2, where region 1 (2) 
is not coupled to the R (L) part of the system via direct 



hoppings or overlaps. Furthermore, regions 1 and 2 are 
connected to each other by ti2 = i?i2 — ESi2- Then we 
may write 



T^(^;) =Tr[Anri2A22r2i], 



gli) and 



(Bl) 



[9ii 



^cc)2i- Similar 



where An = i{gii 
[ESii — Hii — (S£)j^j^] ^ are the spectral density and 
the Green's functions of region 1 in the absence of ti2, 
2^12 = ti2 + ti2G'2iti2, and G21 
expressions hold for A22 and T21 . 

In our case, the regions 1 (2) are made up of all 
atoms in the first (second) phenyl ring plus the sulfur 
and three gold atoms to the left (right) in region C 
(Fig. [2]). To simplify the discussion, we consider the 
electronic structure of the molecule in the junctions as 
separable into a and tt valence electrons, a procedure 
called TT-electron approximation^^ Furthermore, we con- 
centrate on the couplings between those 2p orbitals on the 
ring-connecting carbon atoms, which contribute to the tt- 
electron system. These are oriented perpendicular to the 
respective phenyl rings, and are thus rotated by the angle 
Lp with respect to each other. The indices 1 and 2 then 
refer to these 2p orbitals, and the matrices in Eq. IjBip 
become scalars. Within an extended Hiickel model -ffi2 
is proportional to the overlap S12 (Refs. l4^l47H48h and 
the scalar coupling element ii2(v) — Hi2{tp) — ESui^p) 
at tilt angle (p is seen to be proportional to cos p. 

Because the Fermi energy of gold is located in the 
HOMO-LUMO gap of the organic molecules (Fig. [H, 
G21 can be assumed to be small at Ep. Therefore 
Ti2{^) « tui'p)- Since the (p dependence of An {A22) 
can be expected to be small, the cos^ ip behavior of the 
zero-temperature conductance follows from Eqs. ^ and 

m 

G^{T = OK) = GoT^{Ep) « |ti2(¥')|' ^11^22. 

Here all energy-dependent quantities are evaluated at 
Ep. 

Small deviations from the cos^ Lp dependence of 
G^p{T = OK) can be expected due to higher-order terms 
in the expansion of T12 or other than the tt-tt couplings 
in ti2- These include for example a-a coupHngs of the 
ring-connecting carbon atoms and next-nearest-neighbor 
couplings between regions 1 and 2. 
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